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O ■ Abstract 

The production of and Z° bosons in collisions of longitudinally polarized 
protons allows to access various combinations of the polarized quark distributions 
from the study of double and single spin asymmetries. We compute the 0(a s ) 
corrections to these asymmetries as function of the vector boson rapidity y in the 
hadron-hadron centre-of-mass frame. Detailed studies of the numerical impact of 
the next-to-leading order corrections are presented. The theoretical uncertainties of 
an extraction of polarized quark distributions from future vector boson production 
data are investigated in detail. 
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1 Introduction 



The experimental knowledge on the spin structure of the nucleon has up to now been 
largely restricted to measurements of the polarized structure function g{' n ' d {x^Q 2 ). 
This structure function probes a particular charge-weighted combination of the polarized 
quark distributions. Its mere knowledge is therefore insufficient to disentangle the contri- 
butions of valence and sea quarks to the nucleon spin and for a further decomposition of 
the light quark sea into different flavours. A recent study of asymmetries in semi-inclusive 
hadron production by the SMC collaboration yielded some additional constraints on 
the valence and sea quark polarization. These data are however still insufficient for a 
simultaneous determination of all different polarized quark and antiquark distributions. 

One of the major future projects in the study of the nucleon spin structure is the 
operation of the Relativistic Heavy Ion Collider (RHIC) at BNL with two colliding po- 
larized proton beams at centre-of-mass energies of y/S = 200 . . . 500 GeV [|. The RHIC 
spin programme is expected to start operation three years from now and will offer new 
information on the polarized parton distributions from a variety of observables accessible 
in hadron-hadron collisions ||]. A recent overview of the physics prospects of the RHIC 
spin programme has been presented in ||, where estimates of the expected asymmetries 
and anticipated errors for the most important processes can be found. At RHIC it will be 
possible to access the polarized quark distributions from the study of double and single 
spin asymmetries in the production of V = W^, Z° vector bosons. 

Up to very recently, these asymmetries have only been studied ||, |], || [jj at lowest 
order in perturbation theory, although it is well known from the unpolarized case that 
perturbative corrections to vector boson production cross sections at hadron colliders are 
sizable P, ^, |H], [11], [12], [13], Hf . The calculation of the next-to-leading order (0(a s )) 
corrections to single and double spin asymmetries as function of the vector boson rapidity 
y in the hadron-hadron centre-of-mass frame and the study of their numerical impact is 
the aim of the present paper. The perturbative corrections to double spin asymmetries 
in vector boson production are identical to the corrections in the longitudinally polarized 
Drell-Yan process and can therefore be easily obtained from [[L5|, [16|] . Corrections to the 
corresponding single spin asymmetries have been first considered in |L7] in the context 
of a soft gluon resummation. We shall present a rederivation of these corrections, which 
follows closely our earlier calculation of the QCD corrections to the double spin asymmetry 
in the Drell-Yan process [IB|. 

The QCD corrections to the total asymmetries (integrated over the vector boson ra- 
pidity y) in Drell-Yan process and vector boson production at RHIC have recently been 
studied in QT8J, |19j. Comparison with these provides a strong cross-check of our results. 

A fully consistent numerical study of spin asymmetries in vector boson production 
at next-to-leading order was until now not possible, as the polarized parton distribu- 
tions were only determined at leading accuracy. With the recently calculated polarized 



two-loop splitting functions 
next-to-leading order from fits pll I 



the polarized distributions can now be determined to 
^3|, E4j to polarized structure function data. Having 



a complete calculation of the two-loop splitting functions available, it is now furthermore 
possible to define consistent scheme transformation prescriptions [EOl E3[ for parton dis- 
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tributions, splitting functions and parton level cross sections at next-to-leading order. 
Possible choices of the factorization scheme in polarized and unpolarized Drell-Yan pro- 



cess are elaborated in great detail in fllg, [19j . In our study, we shall only work in the MS 



scheme (with an appropriate correction for spurious terms generated due to the chosen 
representation of 75), which is identical to the scheme used in the recent determinations 
of polarized parton distributions at next-to-leading order PTl E2 . 



In the context of the present study, we shall only consider the rapidity distributions 
of the produced vector bosons, not of their decay products. In the case of Z° boson 
production, identified from the dilepton decay mode, a reconstruction of the boson rapidity 
is indeed possible, such that a direct comparison of our results with experimental data 
is feasible. W ± bosons can on the other hand only be detected from their lv\ decay 
mode, where just the lepton is observed, thus rendering a direct reconstruction of the 
W ± rapidity unfeasible. The calculation presented here for the W ± bosons is therefore 
only a first step towards a prediction of the lepton plus missing transverse momentum 
distribution in polarized hadronic collisions at next-to-leading order. 

This paper is structured as follows: in section |], we define the single and double 
spin asymmetries accessible in vector boson production at RHIC and review briefly their 
phenomenology. Section ^| contains the calculation of the QCD corrections to these asym- 
metries. As this calculation is technically very similar to the calculation of the corrections 
to the polarized Drell-Yan process presented in our earlier work Jl6[ , we will only point out 
the differences to the Drell-Yan calculation without presenting all details of our deriva- 
tion. The numerical magnitude of the corrections will be studied in detail in section 
Finally, section |5| contains a summary of our results and concluding remarks. 



2 Spin Asymmetries in Vector Boson Production 

The production of vector bosons in hadronic collisions can be described accurately in the 
narrow- width approximation (e.g. as production of an on-shell massive vector boson 
by quark-antiquark annihilation. Like in the case of the polarized Drell-Yan process, one 
can study the double spin asymmetry as a function of the rapidity of the vector boson in 
the hadronic centre-of-mass frame: 

(1) 

+ da-~) , (2) 

+ da") , (3) 

with (+) and (— ) denoting positive and negative hadron helicities. 

The unpolarized cross section in the denominator of this asymmetry can be predicted 
using the known unpolarized parton distribution functions. It is evaluated at next-to- 
leading order for pp collisions at = 500 GeV in Fig. [3], using the unpolarized parton 



where 



dAa LL / da 
dy / dy 



dAa LL = ~ f da ++ - da + ~ - da 
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da = - f da++ + da+~ + da" 
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Figure 1: Unpolarized vector boson production cross sections in proton-proton collisions 
at RHIC (VS = 500 GeV). 



distributions from MRS(set A') [g(J and GRV p7[ . The rather small discrepancy between 



both predictions seems to suggest that the unpolarized cross section can be predicted 
reliably. This is however not true in reality, since the ratio of W + and W~ production 
cross sections in proton-proton collisions is proportional to the ratio of the unpolarized 
light antiquark distributions d(x,Q 2 )/u(x,Q 2 ), which is only poorly known at present. 
In principle, a measurement of the unpolarized W boson production cross sections at 
RHIC could help to determine this ratio, as suggested in 0, However, by the time 
RHIC comes into operation, this ratio should already be well constrained from measure- 
ments of the Drell-Yan lepton pair production in proton-proton and proton-deuterium 
collisions [^9], [3(J, thus enabling a reliable prediction of the unpolarized cross sections. 

The phenomenology of the double spin asymmetries is best understood at leading 
order (and neglecting the CKM mixing factors \Vij\ for W ± production), where they can 
be expressed || as ratios of polarized and unpolarized quark distributions 

w+ . A«(x°, M^)Ad(x° 2 , Ml) + Ad(xl M^)Au(xl Ml) 



4Tl(v) 



u(xl M&)d(xl Ml) + d« Ml)u(x° 2 , Ml) 



w Ad(xl MpAujxj Ml) + Aujxl MpAd(x 2 , Mjj 



d(x1, Ml)u(x° 2 , Ml) + u(xl Ml)d(xl M, 



w J 



o EM + a 2 q ) [Ag(s°, Mf ) Aqjxj M 2 ) + Ag(s°, Mf) Aqjxj M 2 )] 

LL[V) " E q (v 2 + a 2 ) [q(xl M 2 )q(xl M 2 ) + q{x\, M 2 )q{x%, M 2 )] 

the distribtions being probed at 



x 



1,2 



M 2 /Se ±y . (5) 



At RHIC with Vd = 500 GeV, these asymmetries are thus sensitive on x lj2 > 0.05, 
where one expects the polarized quark distributions Au(x, Q 2 ) and Ad(x, Q 2 ) to be dom- 
inated by contributions from valence quarks, which are reasonably well constrained from 
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inclusive and semi-inclusive asymmetries in lepton-hadron scattering. The above double 
spin asymmetries can therefore be used as a direct probe of the polarized antiquark dis- 
tributions. Even the flavour structure of the polarized quark sea can be revealed, since 
W + production probes only Ad, while W~ production is particularly sensitive on Au. 

In addition to the above double spin asymmetries, it is also possible to define non- 
zero single spin asymmetries, which are induced by the parity violating couplings of the 
vector bosons. These asymmetries are absent in the (parity conserving) classical Drell- 
Yan process and their measurement requires only one of the incoming hadrons to be 
polarized. Taking hadron 1 (moving in the +y direction) to be polarized and hadron 2 to 
be unpolarized, one can define the following single spin asymmetry 

where 

dAa L = - (da ++ + da + ~ - da~ + - da—) (7) 

and da as in eq.@ above. 

At leading order (again neglecting the CKM mixing factors for production), one 
finds the following simple expressions for these single spin asymmetries || |7| : 

- An« Mpd{x 2 , Ml) + Ad{xl Ml)u(x° 2 , Ml) 
u(xl Ml)d(x° 2 , Ml) + d(x° 1} Ml)u{x% Ml) ' 

-Ad(s°, M^)u{x% Ml) + Au[x\, Ml)d[x\, Ml) 
d(xl Ml)u(x° 2 , Ml) + u(xl Ml)d(x° 2 , Ml) ' 

gg(j + aj) [-Ag(x°, M 2 )q(x° 2 , Mf) + Ag« M 2 )g(x° 2 , M 2 )} 

E q (v 2 + a*) [q(a», M 2 )q{x\, Mf ) + q{x\, M 2 )q{x\, Mf )] " [ ] 

At large and positive rapidity y, (x® > x 2 ), the above expressions are dominated 
by the first term in the numerator, since Aq(x®,Q 2 ) ^> Aq(x\,Q 2 ) for large x\. The 
second term in the numerator is dominant for large and negative y, (x® < x 2 ), since 
q(x 2 , Q 2 ) ^> q{x 2 ,Q 2 ) for large x 2 . In practice, this means that a measurement of these 
single spin asymmetries at RHIC with = 500 GeV can probe H, [7j the polarized 
quark distributions for x > 0.2 from Ai(y > 0) and the polarized antiquark distributions 
for x < 0.12 from Ai(y < 0). A discrimination of different quark flavours is again possible 
from a study of W + and W~ boson production. 

Using the double and single spin asymmetries A^^y), A^y), A^—y) together with 
the unpolarized cross section da(y)/dy, one can calculate the cross section dcr AlA2 (|/)/d|/ 
for any combination A1A2 of hadron helicities. It is therefore possible to express the 
parity violating double spin asymmetries A^iy) and A LL (y) discussed in by linear 
combinations of the asymmetries defined above. 



Af{y) 
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3 Perturbative corrections to the polarized vector 
boson production cross sections 

The derivation of the 0(a s ) corrections to spin asymmetries in hadronic vector boson 
production is in many aspects similar to the calculation of the 0(a s ) corrections to the y 
distribution of lepton pairs produced via the Drell-Yan process in collisions of longitudi- 



nally polarized hadrons. We have recently presented a complete derivation [16] of these 



corrections to the longitudinally polarized Drell-Yan process, which was following closely 



the unpolarized calculation of Altarelli, Ellis and Martinelli [|10 |. 

The only lowest order contribution to the Drell-Yan process comes from the annihi- 
lation of a quark-antiquark pair into a virtual photon 7* (or on-shell vector boson V) of 
invariant mass M 2 . At next-to- leading order, one has two different contributions: the 
0(a s ) correction (emission of a real or virtual gluon) to the qq annihilation process and 
the quark-gluon compton scattering process qg —* <77*(V), where q can be a quark or 
antiquark. Both next-to-leading order parton level cross sections contain divergences as- 
sociated with collinear singularities in the initial state. These divergences can be made 
explicit in dimensional regularization with d = 4 — 2e, they are process independent and 
can be factorized into the bare parton distribution functions. In the case of the Drell-Yan 
process at 0(a s ), all infinities can be absorbed into the bare quark distribution. 

The use to dimensional regularization to compute spin dependent quantities faces 
the problem of representing 75 in n 7^ 4 dimensions. We use the 75-prescription 
of 't Hooft, Veltman, Breitenlohner and Maison (HVBM), which has been consistently 
used in the derivation of the next-to-leading order corrections to the polarized splitting 



functions |2tJ. In this prescription, one restricts the anticommutation property of 75 to the 
physical four dimensions, while 75 commutes in the remaining n— 4 dimensions. The major 
drawback of this formalism in the MS-scheme is the non-conservation of the flavour non- 



singlet axial vector currents due to a non-vanishing first moment of the corresponding 
non-singlet NLO splitting function AP qq ^ + . To restore the conservation of this non-singlet 



axial vector current, a further scheme transformation [33J of the results obtained in the 



HVBM formalism is needed [20]. An explicit formula for this scheme transformation, 
represented by a finite counterterm in the bare polarized quark distribution is given for 
example in [jl6| . 

Inclusive vector boson production in unpolarized hadron collisions is up to 0(a s ) in 
perturbation theory identical to the Drell-Yan process, differences due to the coupling 
of the Z° boson to a closed quark loop || occur only at 0(a 2 ). It should therefore be 
expected that the 0(a s ) corrections to the double polarized cross section dAaiL{y)/dy are 
identical to the corrections to the y distribution of lepton pairs produced in the polarized 



Drell-Yan process [15|, [16|]. This identity could only be spoilt if the vector and axial- vector 
components of the vector boson couplings {v q and a q ) to polarized quarks would acquire 
different higher order corrections due to the breaking of the anticommutation property of 
75 in n 7^ 4 dimensions. 

Using the HVBM representation of 75, one can easily show that it is still possible to 
factor out the combination v q + a q from the squared next-to- leading order matrix elements. 
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The resulting matrix elements then coincide with the matrix elements appearing in the 
polarized Drell-Yan process up to terms which yield only contributions of 0(e) to the 
final result. The 0(a s ) corrections to the double polarized cross section @ are therefore 
indeed identical to the corrections to the polarized Drell-Yan process, as naively expected. 

The single polarized cross section defined in (^) vanishes due to parity conservation 
in the 'classical' Drell-Yan process, which is mediated by a virtual photon. Only the 
parity violating couplings of the and Z° bosons to quarks induce non-zero single spin 
asymmetries. In the lowest order annihilation process, one has to distinguish between the 
annihilation of a quark from the polarized hadron with an antiquark from the unpolarized 
hadron (Aq + q) and the annihilation of an antiquark from the polarized hadron with a 
quark from the unpolarized hadron (Aq + q). The e-dependent parts of these lowest 
order cross sections, relevant for the finite terms generated in the mass factorization, are 
non-identical if the HVBM prescription for 75 is used: 



Aa L (Aq + q -> V) ~ -(1 - e) 



Aa L (Aq + q -> V) ~ (1 + e) . 



These normalization factors of the lowest order cross sections determine the terms that 
have to be factored out from all matrix elements of the 0(a s ) subprocess contributions. 
The factor — (1 — e) is required for the (Aq + q) annihilation process and for the (AG + q) 
and (Aq + G) gluon Compton scattering processes, while (1 + e) must be factored out 
from the (Aq + q) annihilation process and the (AG + q) and (Aq + G) gluon Compton 
scattering processes. Once the correct normalization has been factored out, the derivation 
of the partonic coefficient functions is identical to the derivation in the double polarized 
Drell-Yan process ||16|| . We shall only quote the final results below. 



The normalization factor common to all polarized and unpolarized cross sections is 



and the coupling factors read 



N 



ttG f M^V2 
3S 



(9) 



\Vij\ for W ± 



Cij 



(vf + af)5ij for Z° unpolarized and double polarized , 



ZvidiSij 



for Z° single polarized . 



(10) 



The unpolarized and polarized vector boson production cross sections at next-to-leading 
order can be expressed in a compact analytic form as convolution of parton distributions 
with partonic coefficient functions: 



da 



3— = iVV, 
dy ^ 



x • 



dxi / dx 2 



Dqcj (X\ , X2, X l7 X 2 ) + ^ 



^<2> -^2' 2~ 
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2tt 9q 



M 2 \ 

i> x 2> — ) G(x 1: fi 2 F ) {qj(x 2 ,n 2 F ) +<z i (z 2 ,/4)i 



^2tt 99 



i N J2 c a I dxi f dx 2 

dy tr ^ 



(11) 



x • 



Dgg (Xi, X 2 , x", X 2 ) + —D q q I Xi, X 2 , X^, X 



„0 „.0 



a 



2tt 



[Ag i (xi,/i|)Ag j (x 2 ,/i|) + A(&(xi, /j l 2 F )Aq j (x 2: fi F )j 

A.G(x!,fip) {Aq j (x 2 ,fi 2 F ) + Aq j (x 2 ,li 2 F )} 

+^AD$ ^i,x 2 ,x?,x°,^ {A qi ( Xl , i4) + Aq^xt, AG(x 2 , /j 2 ^, (12) 



o o 

X 2 ,Xi,X 2 , „ 



dAa L 
dy 



i,j Jx l Jx 2 



X • 



Dqij (x±, X 2 , Xi, X 2 ) + n _Dqq ( Xi , X 2 , Xj , X 2 , 2 



a 



s nW 



M 2 



2tt 



X 



{ - Aqi{x u i/ F )qj(x 2 , i/ F ) + A^Xi,/^)?^,/^)} 



o o 

, j/ 2 , j/^, a, 2 5 2 



AG{x u n 2 F ) {gj(x 2 ,/i|) - g.,(x 2) /40} 



2n 99 



ry ry> rf ,-y»0 

xi, x 2 , x 2 



M 



{-Aq^Xi, /j 2 f ) + Aq^xy, fi F )} G(x 2 , fi F ) \ . (13) 



The partonic coefficient functions in the above expressions are identical to the coefficient 
functions occurring in the ^-distributions of the polarized and unpolarized Drell-Yan 
process JTT], [T5], [TT], They read: 

Dqq ■ (Xi, X 2 , X°, X°) = 5(xi — X?) 8{x 2 — X°) , 



(14) 



Dqq ( Xi , X 2 , X-^ , X 2 , ^ ^ 



M 2 " 

/4 



\7l 2 



C F l 8{x x - x?) 5{x 2 -x° 2 ) — - 8 + 2Li 2 (x°) + 2Li 2 {x° 2 ) 



+ ln 2 (l - x°) + ln 2 (l - x°) + 2 In 



In 



1 - X? 1 - xl 



+ [S(x 1 -X? 



x 



2 



, 3/ O 

- In- 



x 2 x 2 x 2 (x 2 — x 2 ) x 2 
x 2 2 + x 2 /ln(l — x 2 /x 2 ) 



~, 2 
A, 2 



\ x 2 — x% 



+ 



/y>0 2 I ^y»2 

X 2 T X 2 



x\ (x 2 - x° 2 )_ 
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Dgq ^Xl, X2, X 1 , X 2 , ^2 



DqJ ixi,X<2,X l ,X 2 , (2 



ADjV [x 1 ,x 2 ,x 1 ,x° 2 ,- 2 - 



where 



M 2 ' 



G A (x 1 ,x 2 ,x° 1 ,x° 2 ) 



H A (x u x 2 ,x\,xl) 



In 



2x° 2 (l-x° 1 ) 
x?(x 2 + x 2 ) 



+ (1 <-> 2) 



I G A {x 1 ,x 2 ,x° 1 ,x 2 y ) A o 



M 2 f 

+ In — <^ S(x! - x\) 5(x 2 - x° 2 ) 



3 + 2 In + 2 In 



XY 



-0 

J/ 2 



= T F 



+ (6(x 1 -x 1 ) 
\8{x 2 - x 2 ) 



, y» 2 I ™ 2 

0\ x 2 ~r x 2 



x\ (x 2 



^ + (1 " 2) )}}- (15 > 



XV 



x° 2 + { Xl - 4j 2 ) In 



0\2\ 



2(xi -x?)(l -x° 2 ) 



+2x° l (x 1 - x?) 



+ 



+ ln 



(x 2 - a%)+ 

M 2 f 5(X 2 - X° 2 ) , n , 2 
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(xi + x1)x 2 
+ ^(xux^x^xl) 



G (xi,x 2 ,Xi,x 2 ) | u -c o (j 



/4 



xf 



Dgg \ X 2 , Xi, X 2 , X 1 , 2 

V A*f 



(16) 
(17) 



T, 



^5(x 2 -x° 2 ) 



X 1 



(2x° -xi)ln 



2(x 1 -x°)(l-x 2 l ) 



+2(xi - x°; 



(x 2 - x°)+ 
, M 2 f <5(x 2 — x2) , n ,11 



AG (xi, x 2 , x°, ^2) 0\ 

- + H (x 1 ,X2,Xi,X 2 ) 
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(18) 
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2(xix 2 + x ( j ) x^)(x ( j ) 2 x^ 2 + x 2 x 2 ) 
x 2 x 2 (xi + x?)(x 2 + x 2 ) 

4x5x 2 {x\x\ + Xix 2 ) 



G c (x u x 2: x 1 ,x° 2 ) 
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// C (xi,x 2 ,x^,x^) 
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The occurrence of double polarized and unpolarized coefficient functions in the single 
polarized vector boson production cross section can be understood by simple arguments. 
The W-boson couples only to one particular helicity of the incoming quarks. The un- 
polarized, single and double polarized qq cross sections receive therefore only a single 
contribution from one particular configuration of quark and antiquark helicities and are 
thus identical up to an overall sign. In the case of quark-gluon Compton scattering, both 
different gluon helicities can contribute; the qg cross sections thus receive contributions 
from two different helicity combinations and the unpolarized and double polarized parton 
cross sections are different. In the case of the single polarized cross section, it is easy 
to identify the polarized quark-unpolarized gluon contribution with the fully unpolarized 
subprocess, and the polarized gluon-unpolarized quark contribution with the fully polar- 
ized subprocess. These identifications are less obvious for the Z° boson production, but 
can be made explicit if the Z° couplings to the different quark handednesses are separated. 

The partonic coefficient functions appearing in the single polarized vector boson cross 
sections have been first studied by Weber in [I7|, where they were obtained as a by-product 
in the calculation of the soft gluon resummation to these single polarized cross sections. 
Our results are in a simpler form than the ones given in |T7| and differ both in the polarized 
gluon-unpolarized quark and the polarized quark-unpolarized gluon subprocesses. The 
discrepancy in the former coefficient can be attributed to a different factorization scheme 



used for the polarized gluon distribution in |Tfl| ; the discrepancy in the latter can only be 
understood to be due to a non-conventional normalization of the number of unpolarized 



gluon states in [17]. A consistency check of our results is given by the quark helicity 
arguments discussed above. These explain the occurrence of the unpolarized and double 
polarized coefficient functions in the expression for the double polarized cross section. 
Moreover, integration of our results over y yields the QCD corrections to the inclusive 
asymmetries, which have been recently calculated by Kamal in [[115]. Expressing the final 
results of [19] in the MS^c-scheme (which is identical to our implementation of the MS 
scheme), we find complete agreement with our results integrated over y. 



4 Numerical results 

The numerical impact of the next-to-leading order corrections derived in the previous 
section can be illustrated by evaluating the asymmetries ([3]) and @ using recent next-to- 
leading order parameterizations of polarized and unpolarized parton distributions. 

When using the polarized GS parton distribution functions [EH], we take = 



231 MeV, the corresponding unpolarized cross sections are then evaluated using the un- 
polarized MRS parton distribution functions set A' f2"6fl . The polarized GRSV distribu- 
tions |22] are consistently used in combination with the unpolarized distributions from 



GRV [^7J and for A„ 4 = 200 MeV. If not stated otherwise, we shall always use /ip = M; 
the strong coupling constant a s is evaluated at fip- All results in this section are obtained 
for longitudinally polarized proton-proton collisions at RHIC (y~S = 500 GeV); the single 
spin asymmetries are always obtained for the configuration where the polarized proton is 
moving in the +y direction. 
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Figure 2: Predictions for double (left) and single (right) spin asymmetries in vector boson 
production in pp collisions at s/s = 500 GeV obtained for different parameterizations of 
the next-to-leading order polarized parton densities. 



It must be kept in mind that all present day parameterizations of the polarized parton 
distributions are fitted only to data on the structure function gi ,d,n (x,Q 2 ). These data 
only provide very loose constraints on the polarized sea quark distributions at large x 
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and are clearly insufficient for a flavour decomposition of the polarized light quark sea. 
Consequently, our present knowledge on the polarized sea quark distributions is very 
poor. This uncertainty is illustrated in Figure |2| showing the vector boson production 
asymmetries A LL (y) and A L (y) evaluated at next-to-leading order using the polarized 
parton distribution functions of GS(A-C) |21] and GRSVs,v [22]. A more extensive study 
at lowest order has recently been presented in ||. 

Based on these present parameterizations, it is not possible to predict sign or magni- 
tude of the double spin asymmetry A LL {y) or of the single spin asymmetry A L (y < 0); 
both asymmetries being dominated by contributions from the polarized sea distribu- 
tions in the nucleon. From the ratio of polarized to unpolarized sea quark distribu- 
tions in present parameterizations, it can only be estimated that \All{%))\ < 0.06 and 
\A L {i) < 0)| < 0.3. Only A L (y > 0) is dominated by the behaviour of the polarized va- 
lence quarks at relatively large x, which is reasonably well constrained by the present 
gi ,d,n (x, Q 2 ) data. Consequently, the different parameterizations yield similar predictions 
for Ai(y > 0), a sizable discrepancy can only be seen for W~ production. This discrep- 
ancy is related to a still sizable difference in the Ad v (x, Q 2 ) distributions in the GRSV 
and GS parametrizations and will most likely be sorted out in the near future when new 
precision data on the deuteron spin structure function gf(x, Q 2 ) will become available. 

The above discussion clearly shows that our present knowledge on the polarized quark 
distributions is still insufficient for precise predictions of asymmetries in vector boson 
production at RHIC. The following studies of the numerical effects of the next-to-leading 
order corrections can therefore only illustrate the qualitative effects of these corrections 
without making definite, quantitative predictions for asymmetries or .fT-factors. 

The contributions of the individual subprocesses (gg-annihilation at leading and next- 
to-leading order and quark-gluon Compton scattering) to the polarized vector boson cross 
sections are shown in Fig. |3|. All polarized subprocess cross sections are obtained with 
the polarized GS(A) parton distributions and are normalized to the full unpolarized cross 
section at next-to-leading order. In the case of the quark-gluon Compton process, we 
distinguish moreover the case where the gluon originates from the proton moving in —y 
direction (qg) and the case where the gluon originates from the proton moving in +y 
direction (gq). This distinction is irrelevant for Aerify), it becomes however important 
for Actx(|/), where the qg process probes Aq(x±, Q 2 ) G(x 2 , Q 2 ) while the gq process is 
sensitive on AG(xi, Q 2 ) q(x2, Q 2 )- 

The numerical impact of the individual contributions to both Aai{y) and A<iLL{y) 
is similar to the impact of these contributions in the ordinary polarized Drell-Yan pro- 
cess [ 16[ . The 0(a s ) correction to the gg-annihilation process enhances significantly the 
lowest order prediction while the quark-gluon Compton process contributes with a sign 
opposite to the annihilation process. A definite prediction of the relative magnitude of 
these two next-to- leading order processes is however only possible for AaL{y > 0), as we 
shall explain below. 

The quark-gluon Compton contributions to Aa^(y) arise mainly from the qg config- 
uration for y > and from the gq configuration for y < 0. Consequently, an extraction 
of the large x behaviour of the polarized quark distributions at next-to-leading order 
from a measurement of Ai{y > 0) will not suffer from the uncertainty on the polarized 
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Figure 3: Contributions of the individual parton level subprocesses to the double (left) 
and single (right) polarized vector boson production cross sections in pp collisions at 
y/s = 500 GeV. All polarized cross sections are normalized to the full unpolarized cross 
section at next-to-leading order. 



gluon distribution, since only the known unpolarized gluon distribution contributes in 
the qg configuration. The numerical magnitude of the next-to-leading order corrections 
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Figure 4: Variation of the double (left) and single (right) spin asymmetries in vector bo- 
son production in pp collisions at y/s = 500 GeV under variation of the mass factorization 
scale. 



to Ai(y > 0) can therefore predicted reliably: inclusion of these corrections enhances 
\Aai(y > 0)| by about 15-25% compared to its lowest order value. 

A quantification of the next-to- leading corrections to A<J LL (y) and Aa L (y < 0) requires 
on the other hand some knowledge on the magnitude of the polarized gluon distribution 
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AG(x, Q 2 ), which enters these cross sections in the quark-gluon Compton scattering pro- 
cess. The same problem was already encountered in the QCD corrections to the longitu- 
dinally polarized Drell-Yan process |16| . It illustrates that any extraction of the polarized 
sea quark distributions at next-to-leading order requires at least some order-of-magnitude 
knowledge on the polarized gluon distribution. 

The change of physical quantities such as cross sections or asymmetries under varia- 
tions of the (unphysical) mass factorization scale \xp enables an estimate of the numer- 
ical importance of unknown higher order corrections. We quantify this uncertainty in 
Fig. f|, showing the double and single spin asymmetries in vector boson production for 
(jLf = 0.5, 1, 2 My. It can be seen that the absolute value of All(v) changes by less than 
0.005 in the central region, the variations become slightly larger in the region where y 
approaches its kinematic limit, but they never exceed 0.01. The variations in Ai(y > 0) 
are even smaller and do not exceed 0.004, they become larger in A^(y < 0), where they 
can amount up to 0.01. All these variations are however significantly smaller than the 
difference between predictions for A LL (y) and A L (y) obtained with various parameteriza- 
tions of the polarized parton distribution functions (cf. Fig. |[). It must be emphasized 
that the relative smallness of these observed variations with the mass factorization scale 
arises partly due to the fact that the scale dependence of polarized and unpolarized cross 
sections largely compensate when taking their ratio in the asymmetry. The unpolarized 
vector boson production cross sections vary by about 15% between fip = 0.5 My and 
Hf = 2M V . 



5 Conclusions 

We have presented a complete calculation of the perturbative 0(a s ) corrections to double 
{All{v)) and single (Ai(y)) spin asymmetries in massive vector boson production in 
collisions of longitudinally polarized hadrons. The results can be expressed in simple 
analytic form and are given in eqs.(|T2|) and (|T3|). These corrections are of particular 
importance for the extraction of the polarized sea quark distributions from measurements 
of vector boson production asymmetries at RHIC. 

The next-to-leading order terms contributing to An(y) agree with the corrections 



to this asymmetry in the longitudinally polarized Drell-Yan process fTo| , |16|| , as naively 
expected. 



Corrections to A^y) had been first considered by Weber |T7[ in the context of a soft 
gluon resummation to single spin asymmetries. While our results for the Aq+q, Aq+q and 
AG + q next-to-leading order subprocesses agree (after suitable change of factorization 
scheme) with ||17|| , we obtain a different result for the Aq + G subprocess (polarized 



quark-unpolarized gluon). This difference could be attributed to a non-conventional 
normalization of the gluonic polarization sum used in [I7[ . A check of our results is given 



by the fact that all single polarized subprocess cross sections can be directly related to 
subprocess cross sections in the double polarized or unpolarized case. This behaviour 
should be expected due to the helicity structure of the vector boson couplings to quarks. 
A further check is the integration of our results over the rapidity y, which reproduces the 
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recent results of [O while showing the same above discrepancy with [IT7J . 



We have demonstrated that the numerical impact of these corrections on the polarized 
vector boson cross sections is very similar to the impact of the next-to-leading order 
corrections in the unpolarized case. The corrections due to the individual subprocesses 
turn out to be sizable. However, a .fT-factor between leading and next-to-leading order 
results can at present only be predicted for A L (y > 0), which is insensitive to the yet 
unknown magnitude of the polarized gluon distribution. 

With the knowledge of the next-to-leading order corrections, it is furthermore possible 
to quantify the uncertainty of the theoretical prediction due to the choice of factorization 
scale and hence the theoretical error on a measurement of the polarized quark distribu- 
tions. We have demonstrated that the absolute value of All(v) varies by less than 0.005 
in the central region under change of the factorization scale. The variations in Ar,(y) are 
only sizable for y < 0, where they amount up to 0.01. These variations are however sig- 
nificantly smaller than the difference between predictions obtained for An{y) and A^(y) 
using different parameterizations for the polarized parton distribution functions. 

To summarize, the calculations presented in this paper allowed us to investigate two 
possible sources of theoretical uncertainty on an extraction of the polarized quark distri- 
butions from double and single spin asymmetries in vector boson production at RHIC. 
We found that the variation of the mass factorization scale had only minor effects on 
all asymmetries, thus indicating that the effect from unknown higher order terms can be 
expected to be small. The effects on these asymmetries from the yet unknown polar- 
ized gluon distribution were on the other hand found to be sizable for both A^^y) and 
A L (y < 0), while being negligible in A L (y > 0). This observation illustrates that a consis- 
tent extraction of the next-to-leading order polarized sea quark distributions from vector 
boson production at RHIC will require at least some order-of-magnitude information on 
the polarized gluon distribution, which is to be obtained from other processes. 
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